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ABSTRACT 

Special  iteration  procedures  for  solving  linear  systems 
of  a  particular  form  are  considered.   It  is  proved  that  the 
iterations  can  always  be  made  convergent  and  some  results  on 
improving  the  rate  of  convergence  are  given. 

The  application  of  these  methods  to  Laplace  and 
biharmonic  difference  equations  is  studied.   In  the  former 
case  it  is  found  that  the  present  schemes  are  inferior  to 
some  previously  proposed  methods.   However,  for  the  biharmonic 
case  they  may  be  of  practical  valuer 
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NYO-8678 

SPECIAL  BLOCK  AND  LINE  ITERATIONS  WITH  APPLICATIONS 
TO  LAPLACE  AND  BI HARMONIC  DIFFERENCE  EQUATIONS 

1,   Introduction.   In  this  paper  we  study  some  iterative 
methods  for  solving  systems  of  linear  equations, 

(1«0)  M  5  =  b  , 

in  which  the  coefficient  matrix  M  is  positive  definite 
and  of  the  special  form: 

M=I-H-V,     HV=VH, 

(la) 

T  T 

H^  =  H  ,  V^  =  V   . 

Systems  of  this  kind  frequently  arise  in  the  finite-differ- 
ence solution  of  elliptic  and  parabolic  partial  differential 
equations.   Two  particular  examples  are  considered  in 
Sections  3  and  l^.  where  the  present  method  is  applied  to  the 
Laplace  and  biharmonic  difference  equations. 

Iterative  solutions  of  (loO)  are  usually  computed  by 
procediires  of  the  form 

(lo2)  N  J^^"^^^  =  P  J^^^  +  b  ,       n  =  0,1,2,...  , 

where      N-P  =  M,       [Nj   f^  0      and     J  is   an  arbitrary  initial 
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guess  at  the  solution.   We  consider  such  procedures  in  which 
the  "splitting"  of  M  is: 

N(a,p,Y)  s  al  -  PH  -  yV  , 
(1.3) 

P(a,p,Y)  =   (a-l)I  -  {p-l)H  -  {y-l)Y   , 

and  the  real  parameters   (a,p,Y)   are  restricted  such  that 

(l.ij.)      a)   |N(a,p,Y)|  7^  0  ,     b)   pY  =  0  . 

The  first  condition  is  of  course  required  to  insiore  that  the 
iterates,   J     ,   are  uniquely  defined  by  (1.2).   The 
second  condition  frequently  insures  that  these  iteration 
equations  are  "easily"  solved,-/   Since  the  matrices  H  or 
V  may  contain  "blocks"  of  elements  of  M  or,  for  partial 
difference  equations,  the  coefficients  of  the  unknowns  on 
various  lines  of  the  mesh  the  present  schemes  include  some 
of  the  so-called  block  [2]  and  line  [Ij  iterations  which 
have  been  studied  previously,. 

It  is  well  known  that  the  iteration  procedure  (1,2) 
will  converge  if  all  eigenvalues,    X    of  N"  P   are  less 
than  one  in  absolute  value..   To  obtain  these  eigenvalues  we 
recall  that,  by  (1,1),   H  and   V  commute  and  are  symmetric 
and  thus  must  have  a  complete  set  of  common  eigenvectors  [3], 
say  \|/j^,   such  that 


-/   If   apY  7^  0   the  solving  of  (1.2)  for  each  iterate  is,  in 
general,  eq\iivalent  to  solving  (loO)  [i»eo   a  =  p  =  y  =  1] 
and  the  iteration  procediore  would  be  ridiculous   in  this 
case, 
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(1.5)  H  iTj^  =  h^  *j^  ,     y\  =  ^^^^' 

Here   h,   and  v,   are  the  corresponding  eigenvalues  of 
H  and  V  and  they  must  be  real  as  a  result  of  the  assxamed 
symmetry.   Prom  the  definitions  (le3-l«>i|-)  we  see  that  the 
\|f,   must  also  be  eigenvectors  of  N   P  and  correspond  to 
the  eigenvalues 


(1.6)     X,  (a,p,Y)  =  1  -  TTT 


pl^-YVj 


Here   p-  =  1-h, -v,   are  the  eigenvalues  of  M,   For  fixed 
(a, P,y)  we  define 


«/   ^   .       I,      ^ 


(1«7)  ^    (a,P,Y)  =  max  |l  -    .  \    , 

k       °-   P^k  '  ^k 

and  the  iterations  (1.2)  will  converge  if  X  (a, P,y)  <  !• 

If  the  scheme  converges  its  rate  of  convergence  [1]  is-^ 

log  l/X   and  thus  the  "best"  scheme  of  the  class  (1<,2-1,1|.) 

is  that  for  which   (a, P,y)  ~  (<i  *P  »Y  )   where 

(1.8)     X**  =  X*(a*,p*,Y*)  =  min  X*(a,p,Y)  » 

a,P,Y 

In  the  next  section  it  is  shown  that  convergent  schemes 
are  easily  found  (Thm«  I)  and  some  results  on  obtaining  a 
best  scheme  are  presented  (Thra.  II )o   Best  schemes,  of  the 
form  (1.2-1,14.),  for  some  Laplace  and  biheoTnonic  difference 
equations  are  obtained  in  Sections  3  and  l\., 

—f    This  definition  of  the  rate  of  convergence  is  valid  only 
if  X*  corresponds  to  a  simple  divisor  of  N~lp.   The 
general  case  is  discussed  in  [1]. 
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2«   Converp:ent  and  Best  Iterative  Schemes.   By  (l.Itb)  we  may 
take  Y  =  0     with  no  loss  In  generality.   Then  with  the 
definitions 


(2.0)  XE-|,     7=  a'     ^^^^)  = 
the  eigenvalues  (l«6)  may  be  written  as 

(2.1)  Xj^(x,y)  =  l-mj^(x)y  . 

For  fixed   (Xj,y)   the  definition  (1.7)  now  becomes 

(2.2)  X*(x,y)  =  max|l-m,  (x)y|  , 

k     ^ 

and  as  in  (1.8)  the  points   (x  ,y  )   are  defined  by 

(2.3)  ^       =  X  (x  ,y  )  =  min  X  (x,y)  . 

(x,y) 

The  range  of   (x,y)   to  be  iised  must  be  determined  such  that 
condition  (l.l^a)  is  satisfied.   This  condition  is  equivalent 
to  the  non-vanishing  of  the  eigenvalues  of  N.   Again  using 
the  common  eigenvectors   ^j^  these  eigenvalues  may  be  fovuid 
and  we  obtain 

(2.1}.)  I  (1-hj^x)  ^   0    for  all  k  , 

as  the  restriction  on   (x,y). 

The  first  result,  which  exhibits  convergent  iterations 
of  the  form  (1.2-1.14-),  can  now  be  stated  as 


THEOREM  I,   For  any  positive  definite   M   of  the  form 
(l»l)  there  exist  open  intervals 

(2.5)  a<x<b,     0<y<  c(x) 

such  that  for  all   (x,y)   in  these  intervals  (2.lj.)  is 
satisfied  and 

X  (x,y)  <  1  . 

Proof.   We  recall  that,  by  the  symmetry  of  H,  V  and 
M  their  eigenvalvies  h.  ,v,  and  ^j_  ,  respectively,  are  real. 
Furthermore,  since  M  is  positive  definite 

(2.6)  p^  =  1-V^k  >  °  » 
and  we  define 

12.7;  h  =  max  h,  ,     iijt  =  min  h,  . 

k   ^        *    k   ^ 

Then  from   (2.6)   and  the   definition   (2.0)      m,  (x)   >   0  ^     for 

all     X     such  that      1-hj^x  >   0,      or  using   (2.7),    for  all     x 

in     a  <  X  <  b     where : 


^*  h 


^^h"*        ^=~¥»  iflih^<05 


(2,a)      S     (b)        a  =  i~  ,        b  =  OG    ,  if     h*  <  0   I 


IT 


(c)        a  =  -aD,      b=—  ,  if     h^>0. 

h 

We   could  also   consider   the    case   in  which     m,  (x)   <  0      in 
some  x-lntervalo      However,    no   such  Interval   exists   if 
h*h^  <  0,      which  seems   to  be   the  most   important   and  more 
general  case. 
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These  limits  are  to  be  used  in  (2«5)  and  then  for  fixed  x 
in  the  indicated  interval  the  eigenvalues  X,  (x,y)   given  by 
(2el)  are  linear  functions  of  y  with  negative  slopes  and 
X,  (x,0)  =  le   Thus  for  y  >  0  we  have   Xj^(x,y)  <  1   and 
introducing 

(2.9)  m  (x)  s  max  m,  (x)  ,    m  (x)  =  min  m,  (x)  , 

k   ^  k   ^ 

we  have 

(2.10)  l-m*(x)y  <  X^U.j)   <   l-m^(x)y  ,    y  >  0  . 

From  the  first  inequality  above   X,  (x^y)  >  -1   if  y  is 
restricted  such  that  l-=m  lx)y  >  =1,   or 

y  <  c(x)  s  ^    . 
Ill  (x) 

Thus  for   (x,y)   in  the  intervals  iZoS),      -1  <  X,  (x,y)  <  1, 
and  (2«l4.)  is  satisfeid  which  concludes  the  proof  of  Theorem  I, 

When  a  value   of     x      in  the   interval    (2,5)   has  been 
determined  a   convergent   scheme   is    obtained  by  taking     y     to 
be   a   sxif ficiently   small   positive   numbero      The   limits   of   the 
x-interval  may  be   estimated  by  means   of   the   usual   techniques 
[I}. J   for  bounding  the  maximum  and  minimum  eigenvalues   of   a 
matrix    (in  this    case      H),      However j    the   rate   of   convergence 
of  a  scheme   determined  in  this  manner  may  be    too   small    (l.ee 
X        too  near  unity)   to  be   of  practical   value,      A   scheme  with 
a  larger  rate   of    convergence  may  be   obtained  by  applying  the 
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COROLLARY.   For  any  fixed  x   In  the  interval  (2.5) 

(2.12)         min   X*(x,y)  =  X*(x,y*(x))  , 
0<y<c(x) 

where 


(2.13)  y*(x)  = 


9 


Proof.   For  x   in  (2,5)  and  y  >  0   the  inequalities 
(2,10)  hold  and  the  equalities  are  each  attained  for  some 
k.   Then  from  the  definition  (2.2) 

X*(x,y)  =  max  i  |l-m*(x)y|  ,  |l-m^(x)y|y  . 

The  two  functions  on  the  right  are  piecewise  linear  in  y 
and,  for  fixed  x,   intersect  at  only  two  points: 

y  =  0  ,     y  =  2[m*(x)+m^(x)j'-'-  >  0  . 

The  minimiom  in  (2.12)  is  clearly  attained  at  the  latter 

point  of  intersection  since   X  (x,y)   is  a  decreasing  ftinction 

of  y  for  positive   y  near  y  =  0.   This  completes  the 

proof  of  the  corollaryc   It  may  be  applied  by  obtaining 

estimates  of  the  maximum  and  minimum  eigenvalues  of  M  and 

then  approximating  m  (x)   and  m  (x). 

If  the  value  of  y  given  in  (2.13)  Is  used  in  (2.1) 

we  have 

^       m  (x)-2mj^(x)+m^(x) 

\(^>7    (x))  =  -^^  .    ,  , • 

m  ( X )  +m^  ( X ) 
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The  definition  (2.2)    then  yields 

m  (x)=m  (x) 

(20II4-)         X*(x)  =  X*(xpy  (x))  =  -^ ^—  , 

m  (x)+ni„{x) 

where  x  is  confined  to  the  interval   (2«5)»   By  the 

corollary  the  best  scheme   (x  ,y  )   defined  in  (2»3)  is 

now  determined^  by  the  value  of  a  single  parameter, 

X  =  X  ,   for  which  (2»li(.)  is  a  minimurao 

To  find  this  best  scheme  we  note  that,  by  the  definitions 

(2oO)  and  (2o9),   m  (x)   and  m^(x)   are  continuous  functions 

of  X   in  a  <  X  <  b   and  they  aj?e  piecewise  differentiable 

(they  are  in  fact  composed  of  arcs  of  hyperbolae).   The 

"corners"  of  these  functions  (ioeo  points  at  which  the 

derivatives  are  discontinuous)  occur  only  at  points  x  =  x 

for  which  m_(x)  =  m  (x)  =  [m  (x)   or  m  (x)]o   Between  any 

two  such  adjacent  corners  we  can  differentiate  (2»li|.)  to 

obtain 

J   3-  o       *ii   m^xy 

(2ol5)        ™  X*(x)  =  ^TTTT-  -^   t~T— J  " 

m   (xj 

If,  in  the  x-interval  in  question^   m  (x)  =  m  (x)   and 

s 

m„(x)  =  m.  (x)   then 
*      z 

,   m  (x)    p.. 
(2,16)        -^  [-1— J  =  T~ 
^  m*(x)    ^ 


W 


Here  and  in  the  remainder  of  this  section  we  neglect 
schemes  for  which  m,  (x)  <  0   (see  footnote  3)«   However, 
it  can  be  shown  that  such  schemes  cannot  be  "besto" 
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Thus  dX   /dx   is  of  one  sign  in  each  interval  between  corners 
Which  implies  that   X  (x)   is  piecewise  monotone  and  takes 
on  its  minimum  value  when  x  =  x   is  at  some  particular 
cornero   The  possible  locations  of  corners  are  easily  fo\ind 
to  be 

^r  s  ^3  r 
If  the  system  (l.O)  is  of  order  n  there  are  at  most 
n(n-l)/2   such  values  of  x   of  which  at  most   (n-l)   of 
them  can  be  corners »   If  the  eigenvalues   p.,   and  h,   are 
known  the  determination  of  x  ^   and  hence  of  the  best 
scheme s,  Is  reduced  to  a  choice  of  one  out  of   (n-l)   or 
fewer  of  the  x    (recalling  that   a  <  x  <  b).   Using  the 

1*3 

definitions  (2o0)  these  results  may  be  summarized  as 

THEOREM  II o   For  any  positive  definite   M  of  the  form 
(1.1)  the  "best"  scheme  (lc2-l,Ii)  with  y  =  0   is  obtained 
for 

(2.18)         a  =  |[m*(x*)-»Tn^(x*)]  ,    p  =  ax*  | 

where     x,      in     a<x     <bj      is    some     x  of   (2el7)   for 

which     X    (jc)     JLs  a  minimum* 
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3.   Laplace  Difference  Equatlonse   The  Laplace  difference 
equations,  in  the  iinlt  square,  have  been  formulated  in  [1] 
for  the  net: 


(3.0) 


X.  =  Iax 


yj  =  JAy 


0  <  i  <  p+1 
0  <  j  <  q+1 


^^  =FT  ^ 


Ay  = 


q+l  • 


In  terms  of  the  unknowns   <|i.    =  4(i-^x,j^y)   and  the  usual 

^»  J 

shift  operators  E    and  E    these  equations  are 

^x         AJ  ^ 

(3.1)   [l-^x(^;i-E.x)-^(ElKy)]  *i.  1  =  ^.,  i  '^~\ 


1  <  J  1  q 


where 


(3.2)   e_  s 


z\y 


'^   2(^x2+^y2)   ' 


©  = 


^x 


y    2(Ax2+..y2)  ' 

^^     ^i,j*  '''o,j  ^  Vl,j  *  '''i^o  ^^  *i,q+l   ^^  assumed 
known.   It  is  shown  in  [Ij  that  the  system  (3.1)  can  be 

written  in  the  form  (l.O)  where!   M  =  I-(L+R)-(B+T)   is  a 

matrix  of  order  pq| 


/-L^+R^ 


(3o3)   (L+R)  = 


0 


\ 


L^-HR^ 


V 


fi 


(B+T)  5 


L  +R  / 
q  q^ 


0  T. 


B^O  T^ 


0 


A 


B,  0 
3 


\0 


.   . q-1 . 

B    0  / 

q 


are   q-order   compound  matrices  whose  non-null  elements   are 
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the   p-order  matrices 


(0 
1   0 


(3oU-)    L^   s  rJ   s 


0\ 


1  ■ , 


0 


1    0 


B.   =  T,    = 


1  <  j  <  q  . 
The  eigenvalues  of  the  matrices  (3,3)  are,  respectively. 


(3.5) 


Pi  -  2®x  °°^  Fi  ' 


r^  =  20   COS  -^  , 


1  <  i  1  p  , 


where  each  p.   has  multiplicity  q   and  each  XT,     has 
multiplicity  p   (see  [Ij  for  proofs). 

To  apply  the  results  of  the  previous  sections  we  take 


H  =  (L+R) 


V  =  (B+T) 


and  conditions  (1,1)  are  then  satisfied.   Let  the  common 
eig^enveators  of  H  and  V  be   \jf.  .,  in  a  convenient  and 
obvious  double  subscript  notation,  with  corresponding 
eigenvalues  h.  ,      and   v.  ..   By  the  multiplicity  of  the 
eigenvalues  (3<>5)  we  have,  as  in  [Ij, 


(3.6) 


^i,j 


Pi  *   "^i,j    ^j  ' 


1  <  i  <  p 

1  1  :!  <  q 


■Lk  - 


which  replace  the   h,   and  v,   of  the  previous  sections. 
The  definition  (2.7)  now  yields,  using  (3.5)  and  (3.6), 

h*  =  p^  ,     h^  =  Pp  =  -  p^  . 

Since   h  h  <  0  ,   (2.8a)  applies  and  the  x-lnterval,  {2,S) , 
of  Theorem  I  is 

(3.7)  -  i  <  X  <  -f  . 

Pi      Pi 

In  this  interval,  with  the  obvious  notation, 

l-p.-r. 

and  with  (3*5)  in  the  above  we  have 

Tl  1  i  1  P 

ra.  t(x)  <  m.  ,(x)  <  m.  „(x)  ,   <  • 

1,1    -  i,J    -  i,q         1  <  j  <  q 

From  these    Inequalities   and  the   definitions    (2.9)    it   is    easily 
foxond  that 

m^   ^(x)    ,      -  — -  <  X  <   T  . -^      , 
■»/     ^    _  <^      P»^  Pi  ~     ^ 


m    (x)    = 


(3.8) 


m„(x)    = 


"^l,q(^^     *       1^1^    <^ 


m,    ,  ( X )    ,      -  — -  <  X  <   T— rr 
1.1  '  D,  —  1-r^ 


Pi  -  1-^1      ' 


P,l  l-t-      -  p^ 


With  these   expressions   in    (2.11)    and    (2«.13)   we    obtain  the 
y-lnterval    (2,$)    of   Theorem   I    and   the    "best"    value    of 
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4t 

y  =  y  (x)   of  the  Corollary, 

However,  the  best  scheme  of  Theorem  II  is  easily  foxmd, 
We  note  that,  by  (3.8),  there  are  only  two  corners  in  the 
interval  (3«7)'   These  occur  at  the  points 


X  =  V  ■■  ■     and    x  = 


l-^^l  1-^1 

Between  these  points  m  (x)  =  m-,      (x)  .   m^(x)  =  m,  -,  ix) 

l,q  *  1,1 

and   forming     X    (x)      as   in   {Zdk.)   yields 

(3.9)  X*(x)    =  t4~  =  ^**   » 

a  constant.   Since  the  best  scheme  is  at  one  of  the  two 
corners,  (3»9)  and  Theorem  II  imply  that  all  schemes  with 

(3.10)  a  =  1—-^  *   P  =  ox  ,   ^rrr^  <  X  < 


are  equivalent  and  hence  in  this  example  there  is  a  one 
parameter  family  of  "best"  schemes.   Since  x  =  1   is 
between  the  corners  we  see  that   a  =  p  =  1  is  in  this 
family  and  yields  a  scheme  which  has  been  called  line- 
Richardson  in  [Ij.   The  principle  eigenvalue  (3<»9)  was 
obtained  there  by  other  means  and  it  was  shown  that  some 
schemes^  not  included  in  the  present  study  are  better  than 
line-Richardson  (for  Laplace's  equation). 

•2/  In  particular  schemes  in  which  all  of  the  improved  values 
CTirrently  available  are  used.   In  line -Richardson  new 
iterates  are  only  used  on  the  line  being  improved  and  the 
old  iterates  are  used  on  the  previous  line<,   The  best 
procedure  has  been  shown  [1]  to  be  one  using  a  combination 
of  old  and  new  values  on  the  previous  line. 
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[].,   Blharmonlc  Difference  Equations,   The  blharmonlc 
difference  equations  on  the  "Interior"  of  the  net  (3.0) 
can  be  formulated  as 

-)2  f^^^^ 

(li.o)  51-e  (e"-^+e,  )-e  (e"-'-+e^  )v  A,  ,  =  b.  .  J 

"+      '      L        X  AX  4x   y  Ay  Ay  j  '''l,j     i,j  )  i<j<q 

Under  certain  boundary  conditions  (given  fxinctlon  and 
second  normal  derivative)  these  equations  with  appropri- 
ately defined  b.  .,   completely  determine  the  numerical 
solution  and  can  be  written  in  the  form  (1«0)  where  now 

(i4..l)  M  =  [I-(L+R)-(B+T)J^  . 

The  matrices   (L+R)   and   (B+T)   are  those  defined  in  the 
previous  section  and  upon  expanding  the  right  hand  side 
in  (i^el)  we  have  a  variety  of  choices  for  the  decomposition 
(1,1),   In  this  section  we  examine  the  selectioni'H/ 

H  =  (L+R)[2I-(L+R)J  , 

ik'2) 

V  =  (B+T)[2I-2(L+R)-(B+T)]  , 


Using  the  common  eigenvectors   'I'-  ^  of   the  previous  section 
the  eigenvalues  of  the  matrices  M,  H  and  V  of  {i\.ol-l].,2) 

—f  With  this  choice  the  equations  (1,2)  which  must  be  solved 
couple  the  unknowns  at  net  points  with  1  =  const*  (i.e. 
horizontal  mesh  lines)  by  means  of  five  term  recursions. 
These  equations  can  be  solved  explicitly  by  a  factoriza- 
tion of  a  5-dlagonal  matrix  into  upper  and  lower  trisingular 
matrices,  in  exact  analogy  to  the  well  known  solution  of 
systems  of  Jacobl ,  or  tri-dlagonal,  form. 
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are   foujid  to  be,    respectively: 

(1J-.3)  h.^    j   =  P^(2-p^)         , 

Applying  the  analysis  of  Section  2  the  definition 
(2,6)  now  yields,  from  (I4..3)  and  {3o5), 

ik.h)  h""  =  h^^j  =  Pi(2-P-l)  ,    h^  =  hp^j  =  -P-l(2+p^)  . 

The  x-lnterval,  (2.5),  of  Theorem  I  is 


(^•^^  p^(iip^)  <  ^  <  Pi(2-p^)   ' 

and  in  this  interval,  from  {k»3)   in  (2,0), 

(l-p.-r;)2    1-p  (2-p.)-v.  . 
i,j     i-p^iz-p^jx    i-p^i2-p^;x 

However,  from  (3*5)  it  may  be  shown  that 


V.    <  V.  .  <  V,  T  , 

and  these  inequalities  in  ([|..6)  imply,  for  x   in   (i|..5), 

/l  <  i  1  P 
m,  , (x)  <  m.  .(x)  <  m.  „(x)  ,  < 
1,1    -  i,j    -  i,q        1  <  j  <  q 

The  definitions  (2,9)  yield,  by  the  above  inequalities. 
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(k,7)    m*(x)  =  max  m.  „(x)  ,   m„(x)  =  min  m.  -.(x)  . 

l<i<p  ^'*1         *      1<1<P   ^'-^ 

It  is  now  not  a  straightforward  matter,  as  it  was  in 
the  last  section,  to  obtain  the  m  (x)   and  m^(x) 
explicitly.   This  is  due  to  the  fact  that  there  are  many 
corners  of  these  functions  (in  fact  p-1   comers  for  each). 
However,  we  can  find  the  best  scheme  of  Theorem  II  without 
a  knowledge  of  the  quantities  that  aj?e  required  to  apply 
the  Corrolary,   To  clarify  the  derivation  of  the  best  scheme 
we  state  the  following  facts  which  are  proved  in  the  Appendix! 

A)  The  corners  of  m  (x)   are  in  the  interval 

— j^ — -r  <  X  <  X   '  where  x     is  the  point  at  which? 

(i^.e8)  m,  „(x)  =  nio  „(x)  . 

For  X   in  the  interval  x    <  x  <  ^  / .-,  ^ — r  , 

-     Pi(2-p^; 

([l.»9)  m  (x)  =  m,   (x)  « 

B)  The  corners  of  m_^{x)   are  in  the  interval 

^(1)^^<^^T^ 

(4.10)  m^  -j^Cx)  =  mg  ]^(x) 

For  X   in  the  interval 


-  j)  -^       '-f 


(i4-.ll)  m^(x)  =  m^  ^(x) 
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C)   All  corners  of  m  (x)   are  to  the  left  of  all 
corners  of  m  (x)  5   ioeo 

We  now  recall,  from  (2«l5"a.l6)  that,  in  the  single 
subscript  notation  of  Section  2g 

ikol3)  sign  [^^^^]  =  sign  [hg-h^J 

where  m  (x)  =  m  (x)   and  m„(x)  =  in.(x).   In  the  present 

S  ■?*•  L 

doxible    subscript   notation   (I4..  13)   becomes,   with  the   aid  of 
(I+.7),   A)^    B)    and   C)s 

sign   i—^]   =   sign/  '    i?'^- 

Prom  the   definitions    (i|.3)    and   {3oB)    the   above   yields 

r<o,     x<x^i^ 

dx*(x)    y  . 

^     (^  >  0  ,     X  >  X(3^) 

and  thus  by  (i^.ol2)  the  best  scheme  must  lie  in  x   <^:^^(i  )<• 
For  X   in  this  interval  we  have,  using  (i4.<,9)  and  (l4.»ll) 
in  (2.11|), 


*/ 


([^..ll^.)   x'^ix)  =  - — i — ^^  =  X*'''  ,   x^^^  <  X  <  X,.. 


Again  X*(x)   is  a  constant  in  the  indicated  interval 
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and  the  best   schemes   form  a  one  psiranieter   family »      Prom 

(i|.olO)  and  (ij.,  13)  in  (2.18)  the  extrapolation  parameters 

corresponding  to  this  family  are 

l-p,(2-p,)+tf  (,) 

(^»15^   ^  =  1-Pl(2-Pi)x    -    P=«^.    ^^^<^<^(1). 

It  should  be  noted  that  the  above  value  of  X    is  a  simple 
function  of  the  corresponding  value,  (3.9),  for  the  Laplace 
scheme j  namely,  in  an  obvious  notationo 

/i  T/\         ,•?»-»■         _      Laplace 
^  Biharmonae   -\+(\  )^ 

Laplace 

Thus  if   0  <  Xt   ,     <  1  we  must  have 
=-  Laplace  — 

"X**  :x'** 

Biharmonic  =  Laplace  ' 


and  the  "best"  biharmonic  scheme  converges  more  slowly  than 
the  "best"  (of  the  class  considered)  Laplace  scheme 0 

Since  better  schemes  than  those  of  section  3  exist  for 
the  Laplace  difference  equations  (see  LI  J)  it  is  reasonable 
to  expect  the  same  to  be  true  of  the  biharmonic  schemes  of 
this  section.   In  fact  the  present  class  of  schemes  do  not 
use  the  latest  improved  value s^   4. .    ,   which  are  known 
off  the  current  mesh-line  which  is  being  improved.   Using 
these  new  values  off-line  might  improve  the  convergence  but 
such  schemes  cannot  be  analyzed  by  the  present  methods, 

I  wish  to  thank  George  H.  Ki±yfor  helpful  suggestions  in 

the  analysis  of  this  section  and  also  for  preparing  the 

Appendix. 
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Appendix 
George  H.  Kuby 

It  has  been  shown  in  Section  l^,    eqs,  ( I4..  3  )- (1+*  7  ) ,  that 

(A.l)     iti*(x)  =  max  m.  ^(x)  ,   m  (x)  =  max  m.  -,(x)  ; 

l<i<p   ^'^        "'      l<i<p   ^'^ 

where   x   Is  confined  to  the  interval 

and  the  required  definitions  are  given  in  (3.5)  and 
(A-3)     m,^^(x)  E  3^  ,  ,,.  H  (i-p,-r.)2  >  0  , 

^ij  =  Pi(2-Pi)  . 

Prom  (A, 3)  we  see  that  each  m. .(x)  is  an  hyperbola  and  thus 
in  the  interval  (A.2),  m  (x)  and  mj^(x)  are  picewise  continuous 
(and  piecewise  monotone)  functions  of  x.   These  functions 
also  have  continuous  derivatives  except  at  certain  points, 
defined  as  corners. 

It  is  clear  from  (A«3)  that 

(A.Ij.)  lim  m-,   (x)  =  +  00  , 

x— >(' TT r  -  0) 

Pl(2-p^ 

and  by  the  continuity  of  m  (x)  and  (A»l)  it  then  follows  that 

(A. 5)  m' (x)  =  ra,  ^(x) 
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in  some  neighborhood  to  the  left  of  the  pole  x   =  l/p,(2-p, ). 
However,  recalling  {3»5)    and  (A, 3)  we  have 

(A. 6)  rn^   (O)  <  m^^-^  q^°^»  1  <  i  1  P-1  » 

and  so  m-,       (x)  must  Intersect  each  ra.   (x),   i  >  1,   at 
i,q  i,q 

some  point  in  0  <  x  <  l/p, (2-p, ).   Thus  m  (x)   ceases  to 

be  m-,       (x)   at  some  point  x     in  this  interval:  this 
l,q 

point  is  the  rightmost  corner  of  m  (x).   The  relation  (A, 5) 
is  valid  in  the  interval  x    <  x  <  l/p-,  (2-p-  ),   which  is 
part  of  statement  A). 

The  point  x    ,   defined  above,  is  the  largest  of  the 
roots  of 

m-i  ^(x)=m.  ^(x),       i  =  2,...,p» 
l,q       l,q  "^ 

In  an  obvious  notation  these  roots  are,  using  (A,3), 

A-p, 
(A»7)  x(p^)  =  Q.p  g  ,  i  =  2,...,p  ; 

where 

(A.8)   A  =  2+2'C'^-p^,   B  S  (2-p^)(l+r-^)^,   C  S  (l+T'^)^-2p£t^^, 

and  by  (3»5)  it  can  be  shown  that 

A>0,    B>0,    C>0,    AC-B  >  0  . 

Consequently: 

^^Pi^   =  AC-B 

dp.       r«.n.  ^2   >  "  ' 


(B-Cp^)' 


-  2k 


and  hence  the   x(p.)   increase  with  p.*   However  the   p. 
increase  with  decreasing  i   and  thus  the  largest   x(p. ) 
Is  x(p^)    =  7i        ,   This  concludes  the  proof  of  statement  A). 
Statement  B)  follows,  in  an  exactly  analogous  manner  and 
we  find  that 


(A. 9) 


r 

A  - 


P2 


(1)  ~   '     • 


where   A  ,  B  ,  C   are  obtained  from  (A, 8)  by  replacing   c  , 

by   "  ^1' 

From   (A. 7)  with      i   =   2      and    (A. 9)   we   form 

/^x       (c-c' )p|-(a'c-ac'+b-b' )p2+(a'b-ab' ) 


X/T   \-x 


^^'  (B    -P2C    )(B-P2C) 

Q(pp) 
(A, 10) 


(B   -p-j^C    )(B-p^C) 

The  numerator  above  is  quadratic  in  pp   and  factors  in  the 
form 

Q(P2)  =  l+r-L[2-p^-P2J[(l-Pi)(l-P2)-r^]  . 

In  this  form  it  is  easily  seen  that   Q(p2)  >  0   (since 

P  2 

T!-,  <  (l-p-i  )    and  p..  >  Pp).   Similarly  it  can  be  shown 

that  the  denominator   in  (A,10)  is  positive  and  so 

^(1)  ^  "^ 
which  concludes  statement  C). 
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